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Abstract. By using the non-relativistic effective Lagrangian approach to bound states, a complete expres- 
sion for the isospin-breaking corrections to the energy levels and the decay widths of kaonic hydrogen is 
obtained up-to-and-including 0(a, rrid — m u ) in QCD. It is demonstrated that, although the leading-order 
corrections at 0{a 1 ^ 2 , (rod — rriu) 1 ^ 2 ) emerging due to the unitarity cusp, are huge, they can be expressed 
solely in terms of the KN S-wave scattering lengths. Consequently, at leading order, it is possible to derive 
parameter-free modified Deser-type relations, which can be used to extract the scattering lengths from the 
hadronic atom data. 

PACS. 11. 10. St - 13.75. Jz - 12.39. Fe 



1 Introduction 

The ongoing DEAR experiment [I] at the DA#NE facil- 
ity (LNF-INFN) aims at an accurate measurement of the 
ground-state strong energy shift and width of kaonic hy- 
drogen, and the strong shift of kaonic deuterium. Prelimi- 
nary results of the measurements for the kaonic hydrogen 
have been reported in Ref. PJ, 



AE{ = 202 ± 45 eV , A = 250 ± 138 eV . 



(1) 



Here, AEf stands for the strong-energy-level shift of the 
ground state of the kaonic hydrogen (total energy shift 
minus certain electromagnetic contributions), and 7\ de- 
notes the width of the ground-state. It should be pointed 
out that these results are in contradiction with the earlier 
measurements BiOJOlKI • see also Fig. [3] below. 

The final goal of the DEAR experiment is to extract 
precise values of the KN S-wave scattering lengths from 
the data by using some counterpart of Deser-type relations 
[7|. Neglecting isospin-breaking corrections altogether, in 
the case of kaonic hydrogen these relations are given by 



AE{-- ri 



2a 3 fj% a,K- 



a K - 



i(ao + ai),(2) 



where fi c denotes the reduced mass of the K~p system, 
and ao,ai stand for the 7 = 0,1 S-wave KN scattering 
lengths in QCD in the isospin limit (a — 0, = m u ). 
In addition, our definition of the isospin limit implies that 
the particle masses in the multiplets are taken equal to 
the charged particle masses in the real world (proton, 7r+, 
K + , • • •). Further, in the experimental proposal it has been 
stated that the precise knowledge of the KN scattering 



lengths could allow one to deduce more accurate values 
for the KN cr-terms and the strangeness content of the 
nucleon. In practice, however, the implementation of the 
above program might pose a rather big challenge to the- 
ory (see, e.g. |H) and references therein). For this reason, 
in this paper we restrict ourselves to the moderate goal 
of relating the KN scattering lengths to the measurable 
characteristics of kaonic hydrogen at the accuracy that 
matches the experimental precision. Using these scattering 
lengths for determining the parameters of the low-energy 
kaon-nucleon interactions is thus out of the scope of the 
present paper. 

It turns out that the isospin-breaking corrections to 
the lowest-order relation given by Eq. are huge. In par- 
ticular, these are much larger than their counterparts in 
pionic hydrogen, or in pionium (typically, a few percent). 
This can be immediately seen, e.g. from the Tabled (see 
below) by comparing the entries in the same column. The 
reason for this qualitative difference will be discussed. In 
addition, the existing predictions in the literature, most 
of which are done in the framework of potential scatter- 
ing theory (see [5lll0lllllll2| for an incomplete list of the 
earlier work on the subject), are anything but consistent 
with each other. In particular, although large corrections 
have been predicted in some of these papers, these effects 
have not been treated systematically - e.g., it is not al- 
ways clear whether all possible large corrections are taken 
into account. Needless to say, all this could make the in- 
terpretation of the results of the accurate measurements 
of the DEAR experiment a difficult task. 

The aim of the present paper is to obtain the for- 
mal relation between the energy shift and the width of 
kaonic hydrogen, and the KN scattering lengths up-to- 
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and-including isospin-breaking effects at 0(a,md — m^) 1 
in QCD, by using the systematic approach to the bound- 
state problem based on the non-relativistic effective La- 
grangians. In the past, this method has already been ap- 
plied to study pionium and kaonium [T31ll4lll5| , as well as 
pionic hydrogen [TBIll7lll8| . It happens to be a very use- 
ful and convenient approach to describe the spectrum and 
decays of this sort of bound states. We shall see that the 
approach is universal: the treatment of kaonic hydrogen 
closely follows the pattern of pionic hydrogen. 



2 Formalism 

As it was already pointed out, the quantities that we 
are aiming to extract from the data on kaonic hydro- 
gen, are the S-wave KN scattering lengths ao, ai eval- 
uated in QCD in the isospin limit, i.e. in the absence of 
the electromagnetic interactions and at m e i — m u . Note 
that we avoid using the threshold scattering amplitude 
calculated with physical hadron masses but in the ab- 
sence of electromagnetic effects, which is sometimes en- 
countered in the literature (see e.g., ["HOT^U] !. The reason 
for this is that this quantity can not be consistently made 
ultraviolet-finite to all orders in Chiral Perturbation The- 
ory (ChPT) 2 . Further, the isospin-breaking effects are pa- 
rameterized in terms of a and — m u . It is convenient to 
introduce a correlated counting of these effects, defining a 
formal parameter 5 ~ a ~ m ( i — m u |16|. The Eq. is 
then valid up-to-and-including 0(5 3 ) in isospin breaking, 
and to all orders in the chiral expansion for the quantity 
ao + eti which is present in this expression. In the present 
paper we modify the relation (J2J), including all terms of 
order <5 7 / 2 and S 4 . 

In order to construct a non-relativistic Lagrangian that 
can describe the spectrum of kaonic hydrogen at 0(6 ), we 
note that: 

i) The only states that are degenerate in mass with the 
K~p state in the isospin limit 5 — > Q, are the states 
K~p + n-f, K°n + n'-f, with n, n! = 0, 1, We explic- 
itly "resolve" only these states in our non-relativistic 
theory, whereas the effect of other intermediate states, 
whose mass is not degenerate with that of the K~p 
state in the isospin limit, is included in the couplings 

1 We use throughout the Landau symbols 0(x) [o(x)] for 
quantities that vanish like x [faster than x] when x tends to 
zero. Furthermore, it is understood that this holds modulo log- 
arithmic terms, i.e. we write also O(x) for xlnx. 

2 Consider the loop contributions to the scattering amplitude 
with both charged and neutral particles running inside the dia- 
grams. The divergent parts which are generated by these loops, 
depend on charged and the neutral particle masses. Since the 
mass difference contains the electromagnetic piece proportional 
to e 2 , in order to cancel all divergences one needs, along with 
the "strong" counterterm Lagrangian, the "electromagnetic" 
counterterms as well. The latter was, however, ruled out from 
the beginning. 



of the non-relativistic effective Lagrangian. 3 In par- 
ticular, the SU (3) breaking scale m s — m u counts at 
0(6°) in our approach. As a result, all effects which 
are non-analytic in the parameter 5 — e.g. containing 
\/5 or In 5, should be produced by the loop expan- 
sion in the non-relativistic theory. To the contrary, the 
generic couplings gi of the non-relativistic Lagrangian 
are regular functions of 6 and can be expanded in Tay- 
lor series 

9i = 9l 0) + agV + (m d - m u ) 5 f + 0{8 2 ) . (3) 

ii) The couplings di that describe the KN scattering in 
the tree approximation (see Eq. (@J) are complex. The 
imaginary parts of the d, can be related through the 
unitarity condition to the transition cross sections of 
the KN initial state into the different inelastic chan- 
nels. In this case, there exist open strong channels - 
e.g. nS, tt°A, etc. The mass gap between these shielded 
two-particle states and the KN state is determined by 
the SU(3) breaking scale. Consequently, the couplings 
di are complex already at 0(6°). This is different from 
the case of the pionic hydrogen, where the imaginary 
part of the effective 2-pion-2-nucleon couplings is of 
order 5. 

iii) The leading strong decay channel in the case of pionic 
hydrogen is ir°n. The phase space for this decay chan- 
nel is proportional to (m p + M v + — m n — M^o) 1 / 2 and 
is thus suppressed by a factor S 1 ^ 2 . For this reason, 
the ratio of the decay widths into the leading electro- 
magnetic channel 777, and into the Tr°n channel counts 
as order 5 1 / 2 only. Numerically, the branching ratio 
into the 777 channel amounts up to ~ 40 % in the 
total decay width. In contrast to this, in the case of 
kaonic hydrogen this branching ratio counts as O(S). 
The measured branching ratio into the leading Ay, 
channels is much less than 1 % [2] (the theoretical de- 
scription of this quantity by using chiral Lagrangians 
[T2Tl2"2"] gives the result consistent with the experiment 
by order of magnitude). Consequently, the perturba- 
tive treatment of the effects due to these channels, as 
it is carried out in this paper, is justified. 4 

iv) Within our approach, it is sufficient to deal with the 
sub-threshold A (1405) resonance indirectly, through 
the (large) KN scattering lengths. The reason for this 
simplification is that the mass gap counts at O(S ) 
in our counting of the isospin-breaking effects, i.e. the 
effect occurs at a "hard scale" and should be included 
in the effective couplings. 

v) The K°n system is a bit heavier than K~p. This simple 
fact has dramatic consequences on the size of isospin- 
breaking corrections, which is nothing but the well- 
known cusp effect (note that the cusp effect is also the 
dominant isospin-breaking effect in some other low- 
energy processes, e.g. in neutral pion photo-production 

3 For the treatment of the 717 intermediate state in pionic 
hydrogen, which is similar in spirit to the approach used here, 
see (Tg|. 

4 We thank E. Oset for interesting discussions on this issue. 
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off nucleons Namely, the loop with the K°n in- 

termediate state at threshold in the non-relativistic 
theory is proportional to [m n +M R o — m p —M K +) 1 / 2 ~ 
\/8 and is real. Its counterpart for the pionic hydrogen 
case is purely imaginary. This means that, in the case 
of pionic hydrogen, the corrections to the real quan- 
tities - the energy shift and width - can not contain 
a contribution from a single neutral loop. Only the 
product of two loops, which is a real quantity, can 
contribute - therefore, the corrections start at order 
5 with respect to the leading-order term. To the con- 
trary, the corrections for kaonic hydrogen can contain 
a single neutral loop. Due to this, the isospin-breaking 
corrections to the Deser formula for kaonic hydrogen 
start at 0{^f5) and are much larger than their coun- 
terparts for pionic hydrogen. 

Despite the differences between pionic and kaonic hy- 
drogen that were discussed above, one may apply exactly 
the same formal approach in both cases to calculate the 
bound state spectra. Below, we closely follow the path out- 
lined in Ref. ^Hi ■ The effective non-relativistic Lagrangian 
is given by 



V 2 

2m T 



V 4 
8mS 



F ecrB D e(X>E - E£>) 



p 2m p p 



s iea(V xE-ExP) 



5m; 



+ J2(K ± )UtD t - M K+ 



D 2 



_ R e(DE-ED) 

+ C K 



+ {K° f \ id t - M R0 + 



2M 
V 2 



K + 



8M^ 



2M R0 8M| 
+ d x ^<ip (K-)^K- + d 2 (^ x (K-yR° + h.c.) 



■■■\K a 



(4) 



Here, F pv stands for the electromagnetic field strength 
tensor (we work in the Coulomb gauge). Further, ip, x, 
and K° denote the non-relativistic field operators 
for the proton, neutron, charged and neutral kaon fields, 
and T>ti> = <9tV> - ieA tp, X»V = V-0 + ieAt/j, DfK* = 
dtK^ieA^K^ , Di^ ± = VK ± ±ieAK ± are the covariant 
derivatives acting on the proton and charged pion fields, 
respectively. The ellipsis stand for the higher-dimensional 
operators and the UV counterterms. The values of the 
couplings Ci can be read off from the matching condition 
for the kaon and nucleon electromagnetic form-factors: 



Fig. 1. Modification of the timelike component of the free 
photon propagator due to the electron vacuum polarization 



M _ 

'K ~ 



M 



Hp, c, 
K+K) 



1 + 2fi p 



$(4 



1 + 2fi f 



where /i p denotes the anomalous mag- 



netic moment of the proton, and 



K 



stand for the 



squared charge radii of the proton and the charged kaon, 
respectively. 

Formally, the vacuum polarization contributions to the 
energy shift and width start at 0(a 5 ). It is, however, well 
known that the electronic vacuum polarization contribu- 
tion is amplified by powers of the factor ^ c /m e ~ 10 3 , 
with m e being the electron mass. It is convenient to count 
the quantity rj = afj, c /m e as 0(6°). Then, the leading- 
order contribution starts already at 0(S 3 ), and the next- 
to-leading order contribution (interference of strong inter- 
actions with the vacuum polarization) comes at 0(6 ). In 
this paper, we explicitly include these contributions in the 
expression of the energy. Note that at the accuracy con- 
sidered here, the whole effect of the vacuum polarization 
reduces to the modification of the timelike component of 
the free photon propagator by an electron loop (see Fig. 1 ) 



A>o(fc) = - 



k 2 



ds 



k 2 3-7T J4 m 2 s + k 



1 / 2m 2 
x - I 1 



1 - 



4m? 



(5) 



3 Energy shift and width of kaonic hydrogen 

The energy spectrum of kaonic hydrogen is obtained by 
using Rayleigh-Schrodinger perturbation theory (for de- 
tails, see Ref. rTHj, see also Ref. for the derivation of 
the energy shift for a generic excited state of the irK atom 
by using non-relativistic Lagrangians). Namely, at the first 
step one constructs the full Hamiltonian H = Ho+Hc+V 
from the Lagrangian given by Eq. (T4J- In this expression, 
Ho stands for the free non-relativistic Hamiltonian of the 
KN pair, He denotes the pure Coulomb interaction be- 
tween K~ and the proton, and the rest of the interac- 
tion is included in the operator V which is treated as 
a perturbation. The general solution of the unperturbed 
Schrodinger equation with the Hamiltonian Ho+Hc, that 
corresponds to the quantum-mechanical Coulomb prob- 
lem for the bound system of spin-0 and spin-1/2 particles, 
is characterized by the quantum numbers n — 1,2, 
j = §, §,•••, m = -j, ■ ■ ■ j and I = j ± ±, 



(Ho + H c )|#„y m (P)) = £„(P)|<F„y m (P)) 
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E n (P) 



M 



K+ 



2(m p + M K +) 2n 2 



= E n + 



2(m p + M K +) 



l^m(P)) 



E 



3 {jm\l{m- s)-s) x 



( 27r )3 v-i-v- -, 2 

**i(m-«)(q)'MM) |P,q, s) , 
|P, q, s) = ft^/iiP + q, s)a f (M2P - q)|0) . 



(6) 



Here, ^ = m p /(m p + M K +), fJ- 2 = M K + /(m p + M K + ) , a 1 *, 
6^ stand for the creation operators for the non-relativistic 
proton and K~, Yi m and (jm\l(m— s)^ s) are the spherical 
harmonics and the pertinent Clebsch-Gordan coefficients, 
respectively, and "P„i(|q|) denotes the radial Coulomb wave 
function, which depends on the magnitude of the rela- 
tive momentum |q|. We further define the "pole-removed" 
Coulomb Green function G n ij(z) and the elastic transition 
operator M„.y(,z) 



P 4 /8Mk* 
--X-- 

+ 



> X > 

p 4 /8nip 



\ / 
\/ 

> I > 

e 



+ 



F D S 
C p , Cp , Cp 



> ii-l > 

% d, 
f 







\ / 
\ .... / 



(I2 d? 



G («) 



1 



-Ho 

G„ y (z) = G(z)- 



G(z) 



1 



E 



z Ho H c ' 

d 3 P \# n l jm (P))(* n l jm (P)\ 



(2tt) s 



^n(P) 



M„y(z) = V + VG nZ j(z)M n y(z) . 



(7) 



By using Feshbach's formalism [25], one can show that 
the (complex) energy shift of the level characterized by 
the quantum numbers njl is given by 

(E n )\& nljm ) + o(S 4 ) , 
(& n i jm {P)\M nl:i {z)\& n i jm (0)) 
= (2tt) 3 5 3 (P) (Vni jm \m nl3 {z)\V nl]m ) . (8) 

In order to evaluate the energy shift up-to-and-inclu- 
ding 0(S 4 ), it suffices to include only the diagrams (a) — (£) 
shown in Fig.|2in the calculation of the operator WL n ij(z). 
The resulting expression is then sandwiched between Co- 
ulomb wave functions. We stress that none of the other di- 
agrams that can be constructed from the Lagrangian (0J , 
neither any possible contribution from higher-dimensional 
operators that are not explicitly displayed in Eq. Q, con- 
tributes to the energy shift at order 8 4 . Here we do not 
provide the details of the calculations, only the final result 
is given. It is convenient to introduce a well-defined split- 
ting of the energy shift into the "electromagnetic" and 
"strong" parts, 5 where the former does not contain the 



5 This naming scheme should not be understood literally. 
For example, the "electromagnetic" contribution depends on 
the electromagnetic radii of the proton and K~ , which are 
determined mainly by strong interactions. On the other hand, 
there are electromagnetic corrections to the couplings d%. 



\ 



/ 

/ 




d 3 d 2 

h 



1 



Fig. 2. The set of diagrams contributing to the energy shift of 
the kaonic hydrogen up-to-and-including 0(S 4 ). Solid, dashed, 
double, dotted, wiggly and spring lines correspond to the pro- 
ton, K~ , neutron, K°, Coulomb and transverse photons, re- 
spectively. The electrons run in the closed loops shown in dia- 
grams (d) and (i). The diagrams (f) and (i) contain Coulomb 
ladders - the contributions with 0, 1, 2, • ■ ■ Coulomb photons 
exchanged. 
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4n 
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7 2 M0?0 725 / ^090 



2tt 



d 2 d 3 2n 



AE^ 



s„(a) = 2(^)(n) - hlna- Inn) , (9) 

n 

and = m n M R a/(m n + Afg-o), go = (2/Uo(m„ + M^o - 
m p - M K+ )) 1 / 2 , V(a;) = r'(x)/r{x). Throughout the pa- 
per, we use dimensional regularization to tame both UV 
and IR divergences. In the above formulae, d stands for 
the number of space-time dimensions and fj, denotes the 
scale of the dimensional regularization. As a check on 
our calculations, we have verified that the electromagnetic 
contributions in the above formulae that come from the 
diagrams (a) + (b) + (c), reproduce the known result of 
Ref. EE] in the limit (r 2 ) = (r 2 K ) = 0. Moreover, with 
Cp = Cp — the above expressions reproduce the energy 
shift in the bound state of two spin-0 particles |24|. 

The contributions due to the vacuum polarization 



AE ( nl ' , AE K n' have been evaluated e.g. in Ref. [27], in the 
same approach as used in the present paper (see also 
|2"5| for the related discussion). At the order of accuracy 
we are working, these contributions do not depend on the 
total angular momentum j. We do not display here the 
(quite voluminous) general result for any n and I. The ex- 
pression for AE^' is given in Eq. (3) of Ref. |27J- Further- 
more, one may write AE^ = -(a 3 /j^/im 3 ) di 5™ c +o(5 i ), 
and relate this quantity to the correction to the bound- 
state wave function at the origin due to the vacuum po- 
larization effect 5™ c = 2 6<Z r n (0)/& n (0). In Ref. |27], this 
correction has been explicitly evaluated for the ground 
state (see Eq. (6) and Table II of this paper), although 
the method used in this paper enables one to make cal- 
culations for any excited level. Here, it is important to 
stress that the "electromagnetic" contributions from di- 
agrams (a) + (b) + (c) + (d), which have to be unam- 
biguously identified and systematically evaluated up-to- 
and- including 0(<5 4 ), are only used for determining the 
so-called "strong shift" (see Eq. J5J) from which the in- 
formation about the KN scattering lengths is extracted. 
Namely, the strong shift is defined as a difference between 
the total energy shift and the electromagnetic shift. In the 
rest of the paper we deal with the strong shift only. 

The equations © do not solve our problem completely: 
the energy shift is expressed in terms of the effective cou- 
plings di which have still to be related to the observable 
quantities. As in Refs. [T51ll6lll8ll^| . this goal is achieved 
by performing the matching for the KN scattering ampli- 
tudes in the vicinity of threshold. In the absence of isospin 
breaking one immediately gets 



j i (o ) = d - 3 (o ) = ^ (ao + ai)) 



(oi-ao),(10) 



where di = d^°' +0(6). However, as one sees from Eq. ||SJ|, 
in order to evaluate AE n ij at 0(6 4 ), d\ should be known 



at 0(6) (for e?2, d% the accuracy of Eq. 1)1 0|l suffices). At 
the required precision, the quantity d\ can be determined 
from matching to the K~p threshold elastic scattering 
amplitude in the presence of electromagnetic and strong 
isospin- breaking effects at 0(6) - the corresponding pro- 
cedure is described in detail in Refs. |lfiU17| . At the first 
step, one removes the one-photon exchange from the spin- 
nonflip part of the relativistic scattering amplitude for the 
process p(p) + K~(q) — > p(p') + K~(q') 



T KN = u(p')\ D(s,t) 



4m r , 



i'^B(s,t)\u(p) 



b l { Sl t) = b(s 1 t)~e 2 F K (t)F 1 (t)( S ~u)/(2m p t), (11) 

where Fx(t) : F\(t) denote the kaon electromagnetic and 
the nucleon Dirac formfactors, respectively, and s,t,u are 
the usual Mandelstam variables. The quantity D'(s,t) is 
singular at threshold, as the magnitude of the relative 3- 
momentum of the proton and kaon in the CM frame |p| 
vanishes. At 0(6), the structure of this singularity is given 
by [TW7] 



e -2 M c (|p|) D'( S ,t) 



IpHo 



Bi 

IpI 
o(IpI) 



Bo In ■ 



T, 



KN 



(12) 



where #c(|p|) denotes the (dimensionally regularized in- 
frared-divergent) Coulomb phase 



MIpI) 



i 



i 



In 



2jp| 



(13) 



In this normalization, the S-wave KN scattering lengths 
and the threshold amplitude Tkn in the isospin limit are 
related by 



T, 



KN 



4tt 1 



+ —)l («o + ai) + 0(VS) . (14) 



The quantity Tkn should be matched to its non-relativis- 
tic counterpart , written in terms of the couplings di. 
A direct calculation with the Lagrangian Eq. which is 
carried out in a similar way as in Ref. |lf>| . yields 



rrNB, 
2 KN 



7 7 2 Mo<7o 
d l - «2 

Z7T 



J 2 7 ( Mo9oY 



- # ^ ( X - 2m) . 



(15) 



The matching condition 2Mk+T^^' = Tkn enables one 
to determine the coupling d\ at the required accuracy. 
Substituting this value of d\ into the expression for the 
strong shift, we finally get the formula in terms of the 



6 



Ulf-G. Meifiner et al.: Spectrum and decays of kaonic hydrogen 



observable quantities, which contains all isospin breaking 
terms up-to-and-including (3(<5 4 ) 



*K - - r n 



2-KM K +n- i 



4?rM 



K + 



x T KN {s n {a) + 2ni) + <5™ 



(16) 



Although Eq. (|16fl formally solves the problem stated in 
the introduction, it is still not well suited for the analysis 
of the experimental data. The reason for this is clear from 
Eq. I|15[l . There it is immediately seen that the unitar- 
ity correction from the K°n bubble (second term in the 
r.h.s. of this equation), whose counterpart in the bound- 
state sector is depicted in Fig. [Sfe, starts to contribute to 
the isospin-breaking part oiT^^ and Tkn at O(VS) (the 
quantity qo is of order V~S). The situation here differs from 
the pionic hydrogen case where the counterpart of the 
quantity go is imaginary because m p + M n + > m n + M^o , 
and the imaginary part of TrN starts at 0(\/~S), not 0(1). 
Consequently, in the ttN case the analog of Eq. (|14f> for the 
real part of the scattering amplitude contains corrections 
at order 5 and not V~5. Exactly the above-described effect, 
which is nothing but the unitarity cusp in the K~p elas- 
tic amplitude, is the source of the huge isospin-breaking 
corrections in the energy shift of kaonic hydrogen, which 
were mentioned in the introduction. 

The above problem can be solved in the following man- 
ner. From Eq. 1)15(1 we see that only the unitarity correc- 
tion to the quantity T^y behaves like y/~S, and all other 
corrections, including corrections in the couplings di, start 
at 0(S) and are regular in S at this order. On the other 
hand, from Eq. l(To")l it is seen that the quantity d 2 can be 
written in terms of the scattering lengths, up to the terms 
of order S. From this one concludes that the corrections at 
O(VS), albeit big, are expressed only in terms of the same 
scattering lengths, that are already present in the Deser- 
type relations at the leading order, 6 and the structure- 
dependent corrections start at 0(5). This counting can 
be implemented in the Deser-formula, e.g. as follows. We 
sum up any number of strong neutral bubbles shown in 
Fig. Ebi since the first term in this expansion contains ex- 
actly the desired singular piece with y/~S. Further, instead 
of Eq. lfTl|) , we write 



T, 



KN 



T 



(o) 



KN 



4^ = ^(1 



2M ( T KN> + 6T KN + 0(6) , 

M K +_ ^ \ (go + a\) + qoa ai 
l + ffi(a +oi) ' 



(17) 



The above equation is nothing but the definition of STkn, 
and our statement amounts to 



ST K n = ati + (m d - m u ) t 2 + o(5) , 



(18) 



6 The same is true for the non-analytic corrections from 
Eq. CHJi which are proportional to a In a. 



where t\,t% are functions of m = | (m u + md), m s and 
Aqcd- For the actual calculation of STkn one has, e.g. 
as in the 7rA^ case, to use the information about the un- 
derlying dynamics which is contained in the low-energy 
effective chiral meson-baryon Lagrangians. With the def- 
inition given in Eq. I|17fl . the formula for the strong shift 
looks similar to that for the case of the pionic hydrogen 



AE' n - l -r n = - , 

"2 2irM K +n 3 



(T K 0) N + 5T KN ) x 



1 - 



afj%s n (a) ^(q) 



AttM 



K+ 



KN 



(19) 



where T K ^, STkn are given by Eq. p7jl. and STkn = 
0(5). This is the final formula, which is best suited for 
the analysis of the experimental data. 



4 Results, discussion and higher order 
corrections 

The isospin-breaking corrections at 0(y/S) that are con- 
tained in the relation of TjPL to the S-wave KN scattering 
lengths ao,ai, are numerically by far the dominant ones. 
These corrections have been derived more than 40 years 
ago 9. 10 by using the X-matrix formalism. However, this 
piece of information should be still supplemented by the 
arguments in favor of the conjecture that the remaining 
corrections are small, i.e. STkn = 0(5), as done in the 
present paper. 

In order to get a feeling how big the corrections to the 
Deser formula can be, we have done a simple exercise. In 
tabled we list the results of the expansion of the quantity 
TjPjL in powers of q a ~ \/S, so that Tn^ denotes T K % cal- 



-(0) 



KN 

= Ti°\ The 



culated up-to-and-including 0(q™) and T KN 
values for the scattering lengths that are needed in these 
calculations, are taken from Ref. 20 .' In these papers, 
the KN scattering amplitudes are obtained by iterating 
tree-level diagrams calculated within ChPT, through the 
Lippmann-Schwinger equation (see also |19ll29j for earlier 
references). In addition, we use the experimental values of 
the scattering lengths given in Ref. [HOj. As we see from 
the table, the corrections at 0(V S) are indeed huge - they 
amount up to a few tens of percent. More precise pre- 
dictions are not possible because the scattering lengths 
themselves are not very well known. On the other hand, 
instead of \ (clq + ai) as contained in the original Deser 
formula, one might determine the combination of these 
scattering lengths T K °^. This combination, which already 

includes the big corrections at 0(v / ^), can be extracted 
from data with a much better accuracy. As one may ob- 
serve from Table ^ the convergence of the expansion in 
the parameter VS is rather good. 

Note also that if one uses the scattering lengths ao = 
-2.24+1.94i, ai = 0.54+0.54i (in fm) as given in Ref. [T3]. 

7 The scattering lengths ao,«i are n °t displayed separately 
in Ref f20) . We thank J. Oiler for providing these values. 



Ulf-G. Meifiner et al.: Spectrum and decays of kaonic hydrogen 



7 





Ref. 20 

a = -1.31 + 1.24* 
ai = 0.26 + 0.66i 


Ref. 30 

a = -1.70 + 0.68* 
oi = 0.37 + 0.60* 


7VTj 0) 
A/"T<°> 
NTf 


-0.52 + 0.95i 
-0.68 + 1.09* 
-0.67 + 1.15* 
-0.65 + 1.16i 


-0.66 + 0.64i 
-0.98 + 0.66i 
-1.04 + 0.73* 
-1.04 + 0.75* 


NT™ 


-0.65 + 1.15* 


-1.03 + 0.76* 



Table 1. Expansion of the KN scattering amplitude in 
powers of qo. Scattering lengths and amplitudes are given in 
fm, and Af = (4tt(1 + M K +/m p ))~ . 



then the convergence of the series in y/S is significantly 
worse. More precisely, one gets: —0.85 + 1.24*, —1.28 + 

1.82*, -1.17 + 2.12*, -1.04 + 2.16*, 1.00 + 2.09* (in 

the notation employed in Table It can be argued that 
this result contradicts the general expectation about the 
size of the isospin-breaking corrections: one sees that in 
these series the corrections at 0(S 3 ^ 2 ) amount approxi- 
mately to 15 — 20 % in the real part of the amplitude. 
Moreover, this result also contradicts our derivation of 
the energy shift: there is no justification for neglecting 
the 0((5 3 / 2 ) terms in the bound-state energy, if in the am- 
plitude their contribution is so large. However, we would 
like to stress that in the approach used in Ref. the 
straightforward introduction of a cutoff to regularize the 
unitarity resummation violates chiral symmetry, since the 
amplitudes are not matched to the chiral expansion 8 . The 
considerably larger values of the scattering lengths than 
those determined by the experiment |30 | , which in its turn 
cause problems concerning the convergence of the series, 
might have been resulted from this implicit violation of 
chiral symmetry. Such a problem does not arise in the 
calculation of Ref. [2]] since an explicit matching to the 
ChPT amplitudes is performed and the regularization is 
done employing by subtracted dispersion relations. That 
this procedure leads to reasonable scattering parameters 
was demonstrated explicitly for the case of pion-nuclcon 
scattering in |32| . 

We further investigate the magnitude of 0(5) correc- 
tions in Eq. I|19fl . The Coulomb corrections that are am- 
plified by In a, are sizable but smaller than those due to 
the unitarity cusp: for the choice of scattering lengths 
from Refs. [20M5Q| , the real part of the correction term 
in the ground state is [9%, 15%], respectively. Again, we 
do not need to know this number very accurately; since 
the Coulomb correction depends on the scattering lengths 
do, a± only, we can use the modified Deser relation which 
includes the Coulomb term to determine 

Below, we briefly consider other corrections contained 
in Eq. |TS)l. The calculation of the quantity STkn proceeds 
analogously to the wN case. This quantity starts at order 
p 2 in ChPT. Further, according to Eq. i|17|) . it is equal 



For a recent discussion on the use of cutoff regularization 
in chiral effective field theories, see |31) . 



to the isospin-breaking part of Tkn at this order. In the 
actual calculations we have used strong (see, e.g. [33113 II 
35 ,36]) and electromagnetic meson-baryon Lagrangians of 
SU(3) x SU(3) ChPT (for the construction principles, see 

E2) 

C 2 = b Q Tr(BB)Tr( X+ ) + b D Tr(B{ X+ , B}) 

+ b F Tr(B[ X+ , B]) + F^dTviQDTviBB) 

+ F§G 2 Tr(Q 2 _)Tr(BB) + F 2 G 3 Tr(BQ 2 + B) 

+ F 2 G 4 Tr(BQ + BQ+) + F 2 G 5 Ti{BBQ 2 + ) 

+ F 2 G 6 Tr(BQ 2 _B) + F 2 G 7 Tr(BQ-BQ_) 

+ F^G 8 Tt(BBQ 2 _) + terms with derivatives , (20) 

where M. = diag(m u , mj, m s ) and Q = 
e diag(2/3, — 1/3, — 1/3) are the quark mass matrix 
and the charge matrix, respectively, B stands for 
the baryon octet field, x+ — 2Bq(uA4u + wAiw), 
Q± = i (uQu^ ± u^Qu), and U = u 2 = exp(-^<£), with 
<P being the pseudoscalar boson octet field. Here, Fq is 
the Goldstone boson decay constant (in the chiral limit) 
and Bq — \(0\qq\0)\/ Fq measures the strength of the 
quark-antiquark condensate (in the chiral limit). Further, 
bo,bD,bf and G,; denote 0(p 2 ) strong and electromagnetic 
low-energy constants (LECs), respectively. Note that 
since our definition of the isospin limit involves the 
physical masses of the charged particles, the terms with 
derivatives of the kaon fields do not contribute to the 
isospin-breaking part of the K~p elastic amplitude at 
0(p 2 ). This quantity is equal to 

st kn = < AA *%u {bo + M + e 2 r _ Gi + G2 _ 2 Gg 

+ 1 -G±- 1 -G 5 + 1 -G (i + \G^, (21) 

where the subscript 'em' denotes the electromagnetic mass 
shift. Using Dashen's theorem, we can replace (Z\M|-) cm 

by (AM 2 ) cnl = AM 2 H at this order. Further, in the 

numerical estimates we use bo = —0.517 GeV -1 , 6_d = 
0.066 GeV^ 1 from Ref. For the unknown electromag- 
netic LECs the order-of-magnitude estimate e 2 Fp \ d\ < 
am p /2-K was used. With these numerical values one gets 
5T KN /T l K % = (-0.5 ± 0.4)% at 0{p 2 ), where T l K ° N = 
M K + 1 Fq is the isospin-symmetric part of the K~p ampli- 
tude at the leading order. The fact that the uncertainty 
in the isospin-breaking part of the KN amplitude turns 
out to be smaller than in the irN amplitude ^7] is not 
surprising: the quantities &Tkn and 5T-kn are of the same 
order of magnitude, whereas the isospin-symmetric part in 
the KN amplitude is increased by a factor M K + /M 7T + and 
additionally by a group-theoretical factor of 2. The results 
are anyway to be taken with a grain of salt: it needs to be 
seen how the results are changed in higher orders in chiral 
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Fig. 3. Predictions of the ground-state strong shift AE{ and 
width Pl . Filled circles correspond to using the Deser formula 
(|2J, empty circles to using T^ instead of | (ao + ai) in this 
formula, and filled boxes to our final formula 119H with STkn = 
8T C = 0. 



The present status of corrections in kaonic hydrogen 
can be summarized by the Eq. (|19(l . Instead of the com- 
bination i (ao + ai) which enters in the original Deser 
formula J2j, we propose to focus on the extraction of 

the quantity T K °^ from the experimental data. The rea- 
son for this is that T K °^ already includes the dominant 
non-analytic corrections in a parameter-free form. The re- 
maining analytic corrections at 0{5) are contained in the 
quantities STkn and <5™ c . The evaluation of STkn within 
ChPT could be interesting, but possibly complicated due 
to the expansion in the strange quark mass. At the present 
stage, in the absence of such calculations, the best is to in- 
clude STkn in the estimate of the systematic error. From 
the above discussion one may hope that the effect from 
STkn should not exceed a few percent, which is a natural 
size of electromagnetic corrections. 



expansion. Finally, we note that the value <5™ c = 0.87% 
given in Ref. suggests that at present accuracy one 
may well include this term in the systematic error. 

In order to visualize the size of the isospin-breaking 
effects, in Fig. [21 we plot the theoretical predictions cor- 
responding to the scattering lengths from Refs. [2T)ll3"n| 
versus the old and new experimental measurements of the 
energy spectrum of kaonic hydrogen PlEBHIHEl • As we 
immediately observe from the plot, the use of the lowest- 
order Deser formula Eq. J2J) can not be justified any more: 
both cusp effect and Coulomb corrections have a size com- 
parable with the present precision of the DEAR experi- 
ment, and should be taken into account during extract- 
ing S-wave KN scattering lengths from the experimental 
data. Note an apparent discrepancy between the recent 
DEAR measurements and the predictions obtained by us- 
ing the scattering lengths from Rcfs. 20,30 . Moreover, 
one sees that when the isospin-breaking corrections are 
applied, the results move away from the DEAR measure- 
ments. We conclude, that the further investigations are 
needed in order to shed light on this interesting issue. 



5 Summary and outlook 

In this paper we derived the formal expression for the 
strong shifts of the energy levels in kaonic hydrogen in 
QCD, up-to-and-including 0(S 4 ) in the isospin breaking 
parameter S ~ a, — m u . The use of the non-relativistic 
effective Lagrangian approach allows one to treat that 
otherwise extremely complicated problem with a surpris- 
ing ease. We discover that large isospin-breaking correc- 
tions arise, in particular, due to the following sources: (a) 
s-channel rescattering with the K°n intermediate state 
(cusp effect), and (b) Coulomb corrections that are non- 
analytic in a. We further prove that the remaining correc- 
tions are analytic in S at O(S). Examining some of these 
corrections, on the other hand, we do not find a big effect 
- the obtained values are at the percent level, which one 
expects to be a typical size of isospin breaking in QCD. 
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